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Abstract. The scenario of homogeneous nucleation is investigated for a first order quark- 
hadron phase transition in a rapidly expanding background of quark gluon plasma. It is 
found that significant supercooling is possible before hadronization begins. This study 
also suggests that spinodal decomposition competes with nucleation and may provide an 
alternative mechanism for phase conversion. 
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1. Introduction 

If quark gluon plasma (QGP) is formed in relativistic heavy ion collisions then 
a phase transition from QGP to normal hadron matter must take place during 
some time of its expansion. Below critical temperature Tc (supercooled), the QGP 
phase becomes metastable and phase transition can be initiated by the nucleation 
of critical size hadron bubbles in a homogeneous QGP phase. The amount of 
supercooling depends on the strength of the transition being more for a strong first 
order transition. As the system supercools, the barrier separating the metastable 
QGP phase from the stable hadron phase decreases and the system may reach the 
point of spinodal instability T5, at which the barrier vanishes. Hence, for a strong 
enough supercooling, the spinodal decomposition may provide an alternate path 
for phase conversion. In case of homogeneous nucleation, the system requires to 
cool upto a temperature Tm at which appreciable nucleation is reached and system 
will be reheated due to the release of latent heat. But if the temperature T,„ is less 
than Ts then the phase transition may proceed via spinodal decomposition. 

For the case of QGP produced in heavy ion collision, the route to hadronization 
either through nucleation or through spinodal decomposition depends sensitively on 
several factors like nucleation rate and expansion scenario. In this work we study 
amount of supercooling and rate of hadronization by solving self consistently the 
nucleation rate along with hydrodynamic equation corresponding to both longitu- 
dinal and spherical expansion. 
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2. Spinodal and supercooling temperature 

As the order as well as the strength of the quark hadron phase transition is still 
an unsettled issue, we consider a more generic form of the potential which covers a 
wide range from very strong to very weak first order phase transition [] 

Vi^,T)=a{T)<j)^~bT^^ + ccb'^, (1) 

where b and c are positive constants which can be expressed in terms of surface 
tension a and correlation length ^. The requirement that at all the temperatures, 
the difference between the two minima should be equal to the pressure difference 
between the two phases fixes the third parameter a{T). The spinodal temperature 
Ts where the quark phase becomes unstable and there exists only one minimum 
corresponding to hadron phase, can be obtained as [] 



Ts 



B 



B + 81a/16i 



1/4 

Tc, (2) 



where B is the bag constant. If the QGP supercools upto this point, it will become 
unstable and may go to hadron phase by spinodal decomposition. For a strong 
enough transition, a/^ is large and Ts is lower as compared to the case when the 
transition is weak. We are interested to know whether the system cools down to 
the temperature Ts- For comparison, we denote the minimum temperature reached 
during the supercooling as Tm- Both Ts and Tm depend on the strength of the 
transition and need to be evaluated properly. While Tg can be estimated directly 
from Eq. (^) , T„i requires a self consistent solution of a set of equations involving 
the nucleation rate and energy momentum conserving hydrodynamical equations 
as described in Ref. [] 

The evolution of the energy density is given by [] 

de ^ Du; _ 47^/3 + C 

dT T 

where D=l for the expansion in (1 + 1) dimension. The factors 77 and C are the shear 
and the bulk viscosity of the medium. For nonviscous plasma (for zero viscosity), 
the above equation follows the Bjorken's scaling solution where T^t is a constant. 
The energy momentum equation needs to be solved numerically for expansion in 
(3 + 1) dimensions. However, retaining the simplicity, we can still use Eq. (^) for 
spherical expansion [], with the choice of D = 3. 



3. Results and Discussions 



The critical temperature is fixed at Tc = 160 MeV. The correlation length is fixed 
at ^ = 0.7 fm and a is treated as a free parameter in the present study. The time tq 
that the plasma takes to cool down to T = Tc will depend on the initial conditions 
and on the expansion dynamics. We assume it of the order of 4 to 8 fm/c. From tc 
onwards, we consider longitudinal expansion if it is 4 fm/c and spherical expansion 
if it is between 6 to 8 fm/c. 
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Spinodal decomposition. 




Figure (|^) shows a plot of Tm/Tc and Ts/Tc as a function of strength of the 
transition, cr/^ for different values of viscosity coefficient rjq at tc=6 fm/c and 8 
fm/c respectively in left and right panels. First consider the case of nondissipative 
plasma with viscosity 77^ = 0. The curves Tm and T5 show a cross over point as 
the transition becomes stronger (large cr/^). For weak enough transition, Tm is 
well above the spinodal temperature Tg. Since depends on expansion rate of 
the medium, the cross over point will sensitively depend on tc] moving towards 
right for the slower expansion (higher tc). If the plasma is viscous [non-zero -qq] 
the Tm is increased further due to slow evolution of the medium. Even the cross 
over point also shifts towards right showing that the nucleation is the dominant 
mechanism over a wide range of a ratios. Figure (§) shows a similar plot for 
a longitudinally expanding system. For comparison the calculation for spherical 
expansion is also shown which gives the lower limit on Tm- One can see that for 
longitudinal expansion, the system remains far from any spinodal instability. 

Similar work has been carried out by Scavenious et al. [], using linear sigma model. 
We use bag model equation of state throughout our work and also incorporate a 
more relevant expansion scenario. 



4. Conclusions 



Although the above study depends on the choice of the parameters, a general ob- 
servation is that for strong enough transition (large cr/^) with zero or very small 
amount of viscosity, Tm < Ts; The system reaches the spinodal instability before 
the amount of nucleated hadron bubbles become significant to begin phase conver- 
sion. The phase conversion in such a case will proceed via spinodal decomposition. 
If the medium is viscous or the transition is weak enough or both, the supercooling 
is much less and Tm > Is; The phase conversion may still proceed through ho- 
mogeneous nucleation. However, depending on the range of the parameters, there 
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ID Expansion (t^ = 4 fm/c) 

o T /T^ for 77 =0 
m' C 'q 




Figure 2. The spinodal temperature Ts/Tc and minimum temperature 
Tm/Tc reached during supercooling 

could be a competition between the homogeneous nucleation and the spinodal de- 
composition as the nucleation and expansion time scales for relativistic heavy ion 
collisions are comparable. 
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